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Abstract 

With the modified Riemann-Liouville fractional derivative, a fractional Tu formula is pre- 
sented to investigate generalized Hamilton structure of fractional soliton equations. The obtained 
qq ■ results can be reduced to the classical Hamilton hierachy of ordinary calculus 
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1 Introduction 

Nobel Laureate Gerardus't Hooft once remarked that discrete space-time is the most radical 
and logical viewpoint of reality. In such discontinuous space-time, fractional calculus plays an 
important role which can accurately describe many nonlinear phenomena in physics, i.e., Brow- 
nian motion, anomalous diffusion, transportation in porous media, chaotic dynamics, physical 
kinetics and quantum mechanics. 

Past decades witness the development of fractional calculus in various fields, such as rheology, 
quantitative biology, electrochemistry, scattering theory, diffusion, transport theory, probability 
potential theory and elasticity. For details, see the monographs of Kilbas et al. [1], Kiryakova 
[2], Lakshmikantham and Vatsala [3], Miller and Ross [4] , and Podlubny [5]. 
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Since Riewe [6] proposed a concept of non-conservation mechanics, fractional conservation 
laws [7], fractional Lie symmetries [8] and fractional Hamilton systems [9-16] have caught much 
attention. In recent study, Fujioka et al found that the propagation of optical solitons can be 
described by an extended nonlinear Shrodinger equation which incorporates fractional derivatives 



Searching for new integrable hierarchies of soliton equations is an important and interesting 
topic in soliton theory. Tu scheme of ordinary calculus [19] is an efficient method to generate 
integrable Hamilton systems. It took various efficient approaches to have obtained many inte- 
grable systems such as AKNS hierarchy, KN hierarchy, Schrodinger system, and so on [20-26]. 
In order to consider the Hamilton structure of fractional soliton equations, some questions may 
naturally arise: (1) Can we have a generalized Tu sheme in fractional case? (2) How to define 
Hamilton equations for fractional soliton hierachy? 

In this study, we start from a Lax pair of fractional order in the sense of the modified 
Riemann-Liouville's derivative [14] and propose a generalized Tu sheme to investigate the Hamil- 
ton structure of fractional soliton evolutionary equations. 

2 Modified Riemann-Liouville derivative 

Generally, there are two kinds of fractional derivatives: local fractional derivatives and nonlocal 
ones. The most used nonlocal operator is the Caputo derivative which requires the defined 
functions should be differentiable. The condition is so strict that many engineering problems 
cannot satisfy, i.e., functions defined on fractal curves, fractional diffusion problem. As a result, 
the Caputo derivative is not suitable for such problems theoretically. 

Several local versions have been proposed: Kolwankar-Gangaks local fractional derivative 
[27-29], Chen's fractal derivative [30, 31], Cresson's derivative [32], Jumrie's modified Riemman- 
Liouville derivative [33] and Parvate's F a derivative [34], among which Jumari's modified R-L 
derivative is defined as 



Here the derivative on the right-hand side is the Riemann-Liouville fractional derivative. The 
nonlinear techniques for such fractional differentiable equations can be found in Refs. [8, 35, 36] 

We can have following results for Jumarie's modified Riemann-Liouville (R-L) derivative. 

(a) The Leibniz product law 

If f(x) is an a order differentiable function in the area of point x, from the Rolle-Kolwankar- 
Jumarie's Taylor series [31], one can have 



[17, 18]. 



r(l - a) J 



< a < 1. 



(1) 



f(x) = lim 



y-tx 



r(l + q)(/(y)-/(x)) 

(y - x) a 



< a < 1. 



(2) 
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Assuming g(x) is an a order differentiable function, the Leibniz product law can hold 



DZ(f(x)g(x)) = g(x)D a J{x) + f(x)D%g(x). (3) 

(b) Integration with respect to (dx) a (Lemma 2.1 of [37]) 

oow ~/V ar^fim = j^—, f /(o«r , o < « < 1. (4) 

(c) Generalized Newton-Leibniz Formulation 

Assume D%f(x) is a integrable function in the interval [a, b}. Obviously, 

f ^—^ J*D2f(x)(dx) a = f(b) - /(a), (5) 

(lT^) Sytfi^T = f(x) - /(a), (6) 



and 

r^jf /(0(«)*-/w. W 

(d) Integration by parts 

With the properties (b) and (c), integration by parts for a order differentiable functions f(x) 
and can be presented as 

— 1— [ b g(x)D«f(x)(dxr = g(x)f(x) \ b a f f(x)D«g(x)(dxr. (8) 

The above properties (a)-(d) can be found in Ref. [33]. 

(e) Fractional variational derivative 

From Jumarie's variational derivative [14] and Almeida' fractional variational approach [15], 
the fractional variational derivative is defined as 

i^ + £(-i>w(^>, ( 9 ) 

where A; is a positive integer. 

(f) From Eq. (2), we can have 

u. m = , im Sll«ig _ Sl^fW, o < „ < i. (io) 

As a result, we can find that 
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3 Fractional Hamilton Structure 



A number of useful attempts have been made to establish fractional variational principles and 
Hamilton system [9-16]. Different types of fractional derivatives may lead to different results, 
for examples, i.e., Baleanu's fractional Hamilton system with Caputo derivative [10], Riemann- 
Liouville type Hamilton mechanics [11], Argwal's Hamilton Formulation with Riesz derivative 
[12] and Jumairie's Lagrange formula [14]. In this section, we revisit Jumarie's fractional Hamil- 
ton system. 

3.1 A Fractional Exterior Differential Approach 

Since Ben Adda proposed the fractional generalization of differential [38, 39], many fractional ex- 
terior differential approaches and applications related to different forms of fractional derivatives 
appeared in open literature [40-42]. A brief review is available in Tarasov's work [43]. 

Starting from the total derivative in the integer dimensional space, and assuming f = f(u,v), 
u = u(x) and v = v(x), where u, v are a order differentiable fuctions and / is a differentiable 
function with respect to u and v, we obtain the total derivative as follows 



On the other hand, if we assume that u, v are differentiable fuctions and / is a a order 
differentiable function with respect to u and v 



df = f u du + f v dv. 



(12) 



Mutiple the both sides of Eq. (12) with T(l + a), we can get 



d a f = f u d a u + f v d a v = f u D«u(dx) a + f v D«v{dx) a . 



(13) 




r(i + a) 




(dv) 



a 



(14) 




or 



2x 2 - a 
r(3-a) 



respectively, if we 



(15) 



3.2 Fractional Hamilton Equations 



We define the fractional functional 




(16) 



Then, we can readily derive the fractional Poincare-Cartan 1-form, which reads 
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u =pd a q- H(dt) a . (17) 

From Eq. (17), we have 

rij-f fib 

d a co = P T ] {dt) a A d a q + d a p A d a q - —d a p A (dt) a - ^-d a q A (dt) a 

op Oq 

azr Pjff 

= [Pt + -HW A d a q + [— (dt) a - d a q] A d a p. (18) 
oq op 

The fractional closed condition d a u> = allows us to obtain the following fractional Hamilton 
equations 



and 



D ? )p= -w (20) 

We must point out, the results Eq. (13), Eq. (15), Eq. (19) and Eq. (20) can be found in 
Ref. [14]. 

4 Fractional Tu Formula and Its Application 

4.1 A Practionalized Tu Formula 

Set A n = A = {a i) j), a it j <G C. Assume A and B € C. Deffine [A, B] = AB - BA. Hence, A n 
is a Lie algebra. The corresponding loop algebra is defined as 

A n = A(n) = AX n , neZ. (21) 

Tu Formula is a beautiful identity to generate integral Hamilton equations. In the past 
decades, many integral Hamilton hierarchies are obtained via this technical [20-26]. Consider 
the fractional compatibility condition, 

<f>W(x,t) = Uct>, 4 \x,t) = Vct>, (22) 

where the fractional derivative is in the sense of the modified R-L derivative [14, 31], and is a 
n-dimensional function vector. 

The compatibility condition of Eq. (22) leades to the generalized zero curvature equation 

_ y{a) + y\ = 0, [U, V] = UV - VU. (23) 
When taking a = (3 = 1, Eq. (22) reduces to the classical zero curvature equation. Set 

n 

U = e (X) + ^2e i (X)u i , {^(A), > i < n} C A n , (24) 
i=i 
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where u = u{u\,u 2 , ...,u n ) T denotes a vector function. By the gradation of B n , define rankA=deg(A), 
then rank (eo(A)) = a, rank (ej(A)) = a%, < i < n are all known. If we take the ranks of Ui as 
a — a-i, \ < i < n, then each term in U is of the same rank a, denoted by rank 

rank (U) = rank (— ) = a. (25) 
If a solution of the stationary zero curvature equation 

- + [U, V] = 0, (26) 

is given by V = J2 m >o HnA _m , (V m )\ = 0, m > 0. rank {V m )\ is assumed to be given so that 
rank {V m )\ = ft, m > 0, then each team in V has the same rank, denoted by 

rank(V) = rank(^) = V . (27) 

Suppose f(A, B) = tr(AB). The follwoing properties can be satisfied 

(a) Symmetry relationship 

f{A, B) = f(B, A); 

(b) The bilinearity can hold 

f(c 1 A 1 +c 2 A 2 , B) = c 1 f(A 1 ,B) + c 2 f(A 2 , B); 

(c) In the sense of the local fractional derivative, the gradient Vsf{A, B) of the functional 
f(A, B) is defined by 

jf(A, B + eC) = f(5 B f(A, B), C),VA, B, C G A n , (28) 

where 5b is variational derivative with respect to B. 

With the fractional variational derivative, we can have the following results, 

S B f(A, B k x a ) = (-l) k 4 ka \ (29) 

where A; is a positive integer and D ka = D^...D". 

v v y 

k 

(d) Communication relationship 

f([A, B], C) = f(A, [B, C]),VA B, C G A n . (30) 
Construct a functional 

W = f(V, U x ) + f(K, - [U, V]), (31) 
where U, V meet Eq. (22), K G A n , rank K=-rank A. 
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With the defined fractional variational derivative, 

5W i \ 5W i \ 

— = V™ - [U, V], W = U X - + [U, V], (32) 
from the above equations, we can derive 

[K, V]^ = [K^, V] + [K, 

= [U X + [U, K], V] + [K, [U, V}} 
= [U x , V] + [[U, K],V] + [[V,U], K] 
= [U x , V] + [U, [K, V}}. 

We can check V' = [K, V] — V\ satisfies Eq. (26) and X- also satisfies Eq. (26) since 
rank(Z)=rank(Vx)=rank(y)-rank(A)=rank(-y ). Therefore, if two solutions of Eq. (23), V and 
V' are linearly dependent, we can get V' = ZV. 

Using Eq. (31) again, we can have a fractional trace identy as follows 

We must point out, the variational derivative here is defined in the sense of the modified 
R-L fractional derivative. 

4.2 Fractional Soliton Hierarchies and Their Hamilton Structures 

Recently, Fujioka et al [17] found the propagation of optical solitons can be described by an 
extended NLS equation which incorporates fractional derivatives. The detailed review can be 
found in Ref. [18]. In view of this point, we consider fractional AKNS hierarchy strating from 
the generalized spectral problem 



^ = U(X,u) = 
Choose a simple subalgebra of A\ 





*= ? " (33) 
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equipped with the commutative relations 

[ei(m), e 2 (n)] = 2e 2 (ra + n) , [ei(m), e 3 (n)] = -2e 3 (m + ra), [e 2 (m), e 3 (n)] = ei(m + n). (35) 
Then, we find that the adjoint representation equation Vx°^ = [U, V] = UV — VU yields 

a o? = 1 c o ~ rb o, bo = 0, c = 0, 

a ix = 1°i ~ rh ii 

= -2b i+1 - 2qa u 
c£) = 2r<n + 2cj+i, i > 1, 

the first few of which reads 

a = -1, 6o = 0, c = 0, 
ai = 0, 6i = g, ci = r, b 2 = ~^i a) > 

1 (a) 1 

We can derive the recurenece relationship 

( c n+1 ) = l(V2- ^q 2rD- x «r ) ( °» ) = L ( °» ) , (36 ) 
V ^n+i / 2 ^ -2gDj«(z -D% + 2qD-« J \ b n J \b n J 

a n = D- a (qc n - rb n ),n = 0, 1, 2, 3 • • • . 

Denoting 

n 

= a(i) ei (n + b(i)e 2 (n - i) + c(i)e 3 (n - i), 
i=0 

Eq. (26) can be written as 

- (V^ n) )i Q) + [U, Vl n) ] = (V M )i^ - [U, V^}. (37) 

It is easy to verify that the terms on the left-hand side in (37) are of degree > 0, while the terms 
on the right-hand side in Eq. (37) are of degree < 0. Thus, we have 

-(V+ ] ) x a) + P, Vi n) ] = 2b(n + l)e 2 (0) - 2c(n + l)e 3 (0). 

Taking an arbitrary modified term for as A n = 0. Notice = it is easy to 

compute, the zero curvature equation 

Uf> - (V^)^ + [U, y (n) ] = 0, (38) 
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which gives rise to 



f,_(qW\ fO -2\fc(n + l)\_ fc(n) 



r 



"' = 1>J, = U oJU+iir^U-) 1- "" 1 '■ (39> 

Here J is Hamiltonian operator. 

For n = 2, we obtain the fractional AKNS equations 



\ Dlr tl = \D«D«r-qr\ 
When a = ft = 1, Eq. (40) can be reduced to the classical AKNS system, 

\ r tl = \r xx - qr 2 . 
In order to use the proposed trace identity, a deirect compute leads to 



(40) 



(41) 



f(V, U x ) = -2a n , f(V, ^) = c n , f(V, ^) = b n . (42) 

X-'xCY (43, 



As a result, 

5u d\ " V b 

Compared with the coefhcients of A~ n_1 



5{ ~ 2an+l) -{l-n)( C :\. (44) 



Su \h n 

Setting n = 1, we can determine 7 = from the initial values. Then we can derive the fractional 
Hamilton function 

_ 2a n+ i (5i?" n _ / c n \ 
n ' <5« \b n J 
The generalized evolutionary equations can be given as 

tn ~ { rfj ) ~ Su ■ (46) 
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5 Conclusion 

Inspired by the previous work [44], in this study we use a different fractional derivative, modified 
Riemann-Liouville derivative, establish a fractionalized Tu sheme for fractional differential equa- 
tions and define a local fractional Hamilton system and derive fractional evolutionary soliton 
hierachies. 

However, there are still other interesting questions needed to be addressed i.e., physical 
mearning of fractional soliton which may be related to fractal media, fractional integral coupling 
method, nonlinear techniques for fractional soliton equations. Such work is under consideration. 
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